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Scanning tunneling spectroscopy applied to the high-Tc cuprates has revealed significant spatial 
inhomogeneity on the nanoscale. Regions on the order of a coherence length in size show variations 
of the magnitude of the superconducting gap of order ±20% or more. An important unresolved 
question is whether or not these variations are also present in the bulk, and how they infiuence su- 
perconducting properties. As many theories and data analyses for high-Tc superconductivity assume 
spatial homogeneity of the gap magnitude, this is a pressing question. We consider the far-infrared 
optical conductivity and evaluate, within an effective medium approximation, what signatures of 
spatial variations in gap magnitude are present in various optical quantities. In addition to the case 
of d-wave superconductivity, relevant to the high-Tc cuprates, we have also considered s-wave gap 
symmetry in order to provide expected signatures of inhomogeneities for superconductors in general. 
While signatures of gap inhomogeneities can be strongly manifested in s-wave superconductors, we 
find that the far-infrared optical conductivity in d-wave is robust against such inhomogeneity. 

PACS numbers: 74.25.Gz,74.72.-h,74.40.-f k 



I. INTRODUCTION 

There is considerable evidence from scanning tun- 
neling microscopy (STM) that some cuprates are in- 
trinsically inhomogeneous on the nanoscale. While 
Bi2Sr2CaCu208+5 (BSCCO) has been extensively 
studiedfii^ii^i^i^i^ii^ spatlal inhomogeneities exist as well 
in La2-a;Sra;Cu04 (Refi^) and in electron doped systems 
such as Pro.88LaCeo.i2Cu04.^'^ On the other hand, the in- 
homogeneities seen in STM can be affected through sam- 
ple preparation mcthods^^. A critical question is whether 
these inhomogeneities exist only on the surface layer to 
which STM is sensitive, or are also present in the bulk. 
Evidence that they are not present in the bulk is provided 
by NMR data on YBa2Cu307_5 (YBCO)^. Also, Loram 
et alji^ have argued that, even in BSCCO, bulk inhomo- 
geneities on the scale indicated in STM experiments are 
inconsistent with specific heat data; although this view 
has been challenged recently by Andersen et alii^ who 
concluded that nanoscale inhomogeneity in BSCCO is a 
bulk property. In view of this conflict, it is clearly impor- 
tant to look for other possible probes of bulk nanoscale 
variation. Such a probe is infrared absorption. In this 
paper, we consider the effect of nanoscale regions on opti- 
cal properties with an aim at identifying their signatures 
in these quantities. 

To understand how nanoscale variations of the super- 
conducting energy gap can affect optical properties, we 
have applied an Effective Medium Approximation (EMA) 
to an inhomogenous system with a distribution of super- 
conducting energy gaps. The merit of this approach is 
that we do not commit ourselves to any particular mi- 



croscopic model of the inhomogeneity, but rather exam- 
ine the macroscopic electrodynamics with a minimum of 
assumptions. In order to consider the case of the high 
Tc superconductors, we include d-wave gap symmetry in 
our calculations and model the gap distribution based on 
BSCCO STM gap maps of Reisi^. As can be seen in 
these and other references, using STM to image a sur- 
face allows for mapping of high and low regions of su- 
perconducting gap magnitude, A. The lattice structure 
of the BSCCO compounds is such that the surface can 
cleave smoothly between adjacent Bi-0 planes, resulting 
in what should be a uniform surface structure. However, 
despite the supposed uniform lattice structure - that is 
to say there is no evidence of restructuring on the surface 
after cleaving - there is a spatially inhomogeneous energy 
gap structure at the surface. These inhomogeneous gap 
structures are commonly displayed in the form of a gap 
map, as seen in Fig. 1(a) In this figure, the blue patches 



are regions with a gap magnitude of about 50 meV and 
the red have gap values of about 20 meV. The largest 
number of regions occur with a gap of ~ 32 meV and 
the variation away from this value is about ±20% as in- 
dicated in Fig. 



1(c) 



The diameters of these regions are 
typically of order 3 nm^i but there is no particular pat- 
tern or underlying structure, rather the regions appear 
to be randomly distributed. It is also important to note 
the discontinuous nature of these gap maps, as can be 
seen in Fig. [21 where a scan along a line across the sur- 
face is shown. Here we see that the gap map contains 
distinct regions with nearly uniform gap values and con- 
sequently does not represent a smoothly varying single 
gap function A(a;), but rather suggests separate patches. 
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FIG. 1: (Color online) (a) STM gap map of surface of over- 
doped (OD) BSCCO with Tc=89K from Ref.2'^(b) Schematic 
of system comprised of spherical particles distributed at ran- 
dom as dictated by histogram (d). (c) Ifistogram of gap dis- 
tribution shown in (a), adapted from Ref.-2-. (d) Re-binning 
of histogram (c) to a set of five gaps. 



each with its own gap value. In implementing the EMA 
as a phenomenological model, we are assuming that these 
STM images directly reflect the bulk gap distribution in 
an attempt to find signatures of this gap inhomogene- 
ity. However, in the following, our results for far-infrared 
optical conductivity will show that no strong signatures 
were observed in mixtures of d-wave superconductors and 
indeed the results of the EMA can be very well approx- 
imated by a single spatially averaged gap value, which 
provides support for the past and continuing use of theo- 
ries and data analyses which assume a spatially homoge- 
neous energy gap in the bulk of high-Tc superconductors. 
It would appear from our work that even if the inhomo- 
geneities are present in the bulk, this assumption will be 
robust. On the other hand, we find that signatures in 
mixtures of s-wave superconductors present themselves 
in the form of a distribution of gap onsets over the range 
of gaps. This distribution results in an optical conduc- 
tivity and optical self-energy which are distinct from a 
single gap system, and that are highly dependent upon 
the range of variation in the gap as well as the actual 
distribution present. In the following section, we begin 
by presenting the theory for the EMA and its use to 
obtain the optical conductivity. In section III, we will 
present our results for d-wave superconductivity based 
on the distribution of gaps seen in STM. Section IV will 
then present results for s-wave superconductivity and we 
summarize our conclusions in section V. 
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FIG. 2: (Color online) Adapted from Ref.^^. Plot of supercon- 
ducting gap value as a function of position along a straight 
line of an STM gap map. 



II. EFFECTIVE MEDIUM THEORY 



It is of interest to us to understand in which systems 
and on what energy scales we might expect to see indi- 
cations of inhomogeneities in bulk optical properties in 
order to make generalizations on, or identify, inhomo- 
geneities in new materials. For this reason, we seek a 
purely phenomenological method of examining these ef- 
fects, and here we will turn to Bruggeman's EM A^^'^^ . 
This EMA has been rigorously tested and previously 
applied to both microscopic metallic and microscopic 
superconducting-normal mixtures^^'i^iii^'^SiSl. It has also 
been used more recently to mode P^i^^ the very low fre- 
quency optical response in the cuprates and to discuss, 
in particular, the THz data on BSSCO from Corson et 
aL— 1^ In certain circumstances, an alternative choice of 
EMA might be used, the Maxwell-Garnet EMA^ how- 
ever, that EMA applies when there are a dilute set of 
regions embedded in a large uniform background, which 
is not the case seen here in the STM of Fig. 1(a) In the 
following, we wish to make a semi-classical approxima- 
tion, treating the system classically as having randomly 
distributed regions, each a uniform superconducting ma- 
terial, with properties defined by non-classical BCS re- 
sults. The implicit assumption here is the randomness of 
the distribution of regions, since it has yet to be under- 
stood how and if these inhomogeneities are distributed 
throughout the bulk material. 



For understanding the optical response of the entire 
system, we now summarize the idea of the EMA. An 
electromagnetic wave, of frequency w, moving through a 
composite material will undergo multiple scattering and 
absorption events. In a bulk material these events are 
accounted for as an average by means of an effective com- 
plex dielectric function em(w). If the bulk material is a 
composite, it is assumed to be composed of a number of 
regions of uniform dielectric functions ei{uj). The scatter- 
ing caused by these regions is dependent upon both the 
shape and size of the regions (grains). These scattering 
events can be taken into account by allowing each grain 
to become polarized, producing its own additional field. 
When we extend this to a macroscopic system, we add 
the constraint that the entire material should have a zero 
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net polarization, i.e 
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m = 0, 



(1) 



where fi is the relative fraction of the occurance of a cer- 
tain type of grain indexed by i and Pi is that constituent's 
polarization. Knowing this, we can solve for the real and 
imaginary parts of em('^) in the following equation: 



(£» - Cm) 

ge^ + {1- g)er. 



0, 



(2) 



where g is the depolarization factor, determined by the 
shape of the constituent grain i^^i^^ The form of the po- 
larization in Eq. ([2]) can be seen as the electric dipolc 
contribution of a multipole expansion. Subsequent terms 
in the polarization (magnetic dipole, electric quadrupole, 
etc.) scale as functions of the grain size. In the spherical 
case these terms scale as {uRi/c)'^ = {Ri/X)^ where Ri 
is the radius of the individual grains and A is the wave- 
length of the field. For this case, the correction to Eq. ([2]) 
to next order would be: 
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We can see certain limits where the second term can be 
neglected. In the small grain limit, we would expect the 
far-infrared (long wavelength limit) response to be com- 
pletely dominated by the first term. As this is the region 
of interest for detecting signatures of the superconduct- 
ing energy gap, the EMA as in Eq. with g set to 
a value of 1/3 (spherical grains), is appropriate for our 
purpose. 

To implement the above EMA, we require the fre- 
quency dependent complex conductivity (t(w) = cri(w) -I- 
i<T2{Lo) as the dielectric function can be written as e{uj) = 
Coo + i[4:Tra{cLj)/uj], where eoo is the high frequency dielec- 
tric constant. The constituent conductivities that we use 
for our purpose here are calculated according to Zimmer- 
mann et al.,^^ for the case of BCS s-wave superconduc- 
tors with variable impurity scattering and according to 
Schiirrer et al.,^^ for the BCS d-wave case. In this paper, 
we consider only an elastic scattering rate in the optical 
conductivity and have not included inelastic scattering, 
as might arise due to electron-boson interactions, such as 
phonons and spin fluctuations. For a system of N differ- 
ent types of superconducting grains, the EMA of Eq. ([2|) 
can be expanded to the form of an A^*'* order complex 
polynomial in em(ijj) or by setting = 1 + i(47r(Tfc/a;) 
{k = 1,2, . . . , N and m) as a complex polynomial in (T.,„. 
For example, in e we now have: 

eZH + aiet\^) + ' • ■ + «ive™('^) - 0, (4) 



where the coefficients vary based on the choice of N . 
For an A^ = 2 system, one can solve directly for e™ by 
substituting in its complex form, tm = eim + i^2m and 
solving for the real and imaginary parts separately. An- 
alytically, this becomes largely unfeasible for N > 2. For 
this reason, considerations of composites, with a large 
number of constituents, N, have not been thoroughly ex- 
plored. Numerical solutions, however, are much more 
feasible. All that is required is solving for the set of a^'s 
for each choice of N. This turns out to be a difficult task 
for A^ > 3 and it is preferable to use a program to solve 
this algebraically. For example, for the N=2 case: 
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7262 (w)(l - g) + /2.gei(w) 



-/iei(^)3e2(t^) - /2e2(w)ffei(w) 
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(6) 



where bo — fi{l ~ g) + /2(1 — .g). We can see that these 
coefficients depend on the dielectric function tiioj) (or 
conductivity) of each constituent, as well as their rela- 
tive volume fractions fi and the depolarization factor, g. 
These are the key pieces of information required to per- 
form such an effective medium calculation, and as such, 
will be explained in detail for each case we wish to ex- 
amine. Knowing the numerical coefficients, in Eq. 31 
reduces the EMA calculation to the problem of numeri- 
cally solving an N*'* order complex polynomial. Such a 
polynomial has N solutions, of which only one is gener- 
ally completely physical as €2(1^) must be strictly positive 
for each choice of lo. 



III. D-WAVE ENERGY GAP 
INHOMOGENEITIES 

In order to determine possible signatures in the bulk 
due to the inhomogeneities seen by STM, we created a 
mixture of superconducting patches which is motivated 
by Fig. 1(a) and shown schematically in Fig. 1(b) We 



tried to maintain the same gap distribution seen in exper- 
iment and shown in the histogram of Fig. 1(c) but have 
re-binned the distribution so that it gives only five frac- 
tions fi in order to reduce the complexity of the numerical 
calcuation. This is shown in Fig. 1(d) Furthermore, we 



assumed that each constituent has d-wave gap symmetry 
similar to what one would expect in the bulk BSCCO 
conductivity^* and have used the BCS d-wave conduc- 
tivity program of Ref.^^ to provide the input conductiv- 
ities to the EMA. The BCS conductivity, while missing 
the inelastic component, does exhibit similar qualitative 
features as seen in experimental results, and can be ap- 
plied with a minimal number of parameters: the energy 
gap. A; the plasma frequency, ujp] the temperature, T; 
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FIG. 3: The real part of the optical conductivity for a d-wave 
superconductor at T = 20K for different gap values chosen to 
match those in the re-binned gap map histogram [Fig, 
The curves shown here have = 2 eV and = 0.2 meV 



and the scattering parameter = 1/{2ttt), where 1/t is 
the elastic scattering rate. Since we are interested in the 
effects of varying only the gap parameter and we wish to 
limit the parameter space, we compute the EMA using 
constituents with similar plasma frequencies and scatter- 
ing rates, but a range of A values. The real part of the 
input conductivities of the constituent superconducting 
grains are shown in Fig. [31 These curves were calcu- 
lated with realistic parameters, as given in Refill, i.e., 
ujp = 2eV ,t+ = 0.2 meV, at T = 20 K. 

The overall characteristics of the resulting EMA in this 
case (shown in Fig. [4]) remain what would be expected of 
a composite of a d-wave superconductor with scattering 
in the clean limit. Indeed, if we superimpose the conduc- 
tivity curve for the average gap value, which dominates 
the gap distribution seen in STM, we find that there is 
not much difference between the two curves and therefore 
we conclude from these calculations that even if the inho- 
mogeneity persists throughout the bulk, the size of vari- 
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FIG. 4: (Color online) Solid black line: EMA calculation in- 
volving constituents described in Fig. |3]with volume fractions 
from Fig. 1(c) (also shown as an inset here). Dashed red 
curve: an overlay of the A = 33.3 meV curve of Fig. [3] for 
comparison. 



ation in the gap, on scales seen in STM, cause very little 
variation in the overall conductivity of the sample from 
that of a single gap system, and therefore it is not pos- 
sible to determine the presence of inhomogeneities from 
bulk far-infrared optical measurements. We note, as well, 
that variation in scattering rates of constituent grains (of 
up to a factor of 10) does not change the qualitative re- 
sults of an EMA mixture with d-wave symmetry. We do 
not expect to have large values of the scattering rate, as 
this would quench the superconducting state, and there 
is no evidence from STM or otherwise for normal regions 
mixed in amongst the regions of inhomogeneous super- 
conductivity at low temperatures. We conclude that the 
narrow range of gap values seen in STM coupled with 
the rather smooth nature of the conductivity for a BCS 
d-wave superconductor makes the observation of possi- 
ble gap inhomogeneities uncertain through far-infrared 
optics. Adding in realistic inelastic scattering would not 
change this conclusion as this type of scattering will not 
introduce new sharp features in the conductivity. As a 
final point, most approaches to calculating optical prop- 
erties of the cuprates assume a homogeneous model. Our 
results indicate that it is appropriate to compare results 
of such calculations for the average gap with experiment, 
and that ignoring inhomogeneities is a valid first approx- 
imation. 



IV. S-WAVE ENERGY GAP 
INHOMOGENEITIES 

With the knowledge of the previous section, we are 
now motivated to find possible systems where gap in- 
homogeneities could have a clear signature in optics and 
indeed, it is clear that a system with a conductivity curve 
which contains sharp features would be an ideal candi- 
date. Consequently, we now examine the effect of a distri- 
bution of gaps on a conventional s-wave BCS supercon- 
ductor. We begin by calculating the zero temperature 
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FIG. 6: (Color online) EM A result for two-grain (left) and 
five-grain (right) composites determined from the histograms 
in Figs. 5(b) and 5(c) The dashed red curve is the single-gap 
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calculation using the average gap value of the distribution. 



FIG. 5: (Color online) (a) Real part of the conductivity for 
different values of the energy gap A, as discussed in the text, 
(b) Histogram of 2 grain mixture, (c) Histogram of 5 grain 
mixture. (1— >5 refer to lowest to highest gap values.) 



BCS s-wave conductivities for two simple distributions 
of gap values shown in Figs. 



5(b) and 5(c) 



We have as- 
sumed a constant elastic scattering rate l/r = 50 meV, 
and plasma frequency LOp — 2eV for all grains and this 
puts the system in the regime where 2 A < 1/t, which is 
a mildly dirty limit. We calculate this EMA mixture un- 
der the assumption that each constituent maintains the 
same overall transition temperature, therefore, variations 
in A value can be seen as analogous to changing the ^ 
ratio from 3.53 7.53 in steps of 1. Curves for the real 
part of the conductivity, used for input to the EMA, are 
shown in Fig. 



5(a) 



In the two grain composite. Fig. [S] (left column), we 
observe that the separation of gap onsets is sufhcient to 
clearly display a two-gap result in the real part of the 
conductivity, ui (top frame). The real part of the con- 
ductivity rises sharply out of zero as it does in all of the 
This is followed by a first peak 



curves in Fig. 5(a) 



af- 
ter which a second rise begins at twice the value of the 
largest gap in our two gap distribution. A second peak 
follows at higher energy before the normal state Drude 
form is recovered. This is fundamentally different from, 
and should not be confused with, a two-band supercon- 
ductor such as MgB2.^^ It should be noted that the EMA 
is distinct in this two-gap system from a fractional aver- 
age of its constituents since both aim and a2m depend 
on both the real and imaginary parts of the constituents. 
In other words, for a two component mixture labeled 
a and &, where (Tq = uia + i<J2a and (if, = an, -f i(72b, 



the EMA results in aim = o-im((Tia, cr2a, ctk,, o-2fc) and 
CT2m = cr2m(CTiQ, cr2a, CTib, cr2h). This additional depen- 
dence of the real part of the conductivity on the imag- 
inary parts of its constituents in general leads to an 
increased low frequency contribution to aim from the 
imaginary conductivities (physically seen as scattering 
from constituents). The imaginary part multiplied by 
cij, namely a;CT2(w), is shown in the lower left panel. As 
the frequency goes to zero, we find the inverse square of 
the London penetration depth for the composite system. 
At finite frequency, there are structures in Loa2{i^) cor- 
responding to the frequencies of the onset of absorption 
of each of the two grains included in the model. These 
structures would be easily identifiable. 

When we consider a less-separated five-grain distribu- 
tion [seen in Fig. 5(c) with results in Fig. [S] (right col- 



umn)] , we immediately recognize the loss of clear onset 
points for each gap value and instead see a broadened rise 
in ai{uj) (Fig. [51 top right frame), which quickly rejoins 
an average gap fit at frequencies just above the highest 
2 A point. Shown here is the EMA calculation as the 
solid black curve with the dashed red curve being the 
calculation for a single gap system using the average gap 
of the distribution as the input gap. It is clear that in 
this more realistic case the signature of inhomogeneities 
is not as readily identifiable as for the two gap case. This 
is also true for the imaginary part, uJa2{^^)^ shown in the 
lower frame; now a single minimum is seen similar to the 
single gap case, but somewhat more broadened. 

The real and imaginary parts of the conductivity are 
not the only optical quantities that enter a more modern 
discussion of the subject. In fact, it has been found to 
be very useful to introduce an optical self-energy Y1°p{uj) 
defined in terms of the generalized Drude formula for the 
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optical conductivity, (y{uj), which is written as 



a{uj) 



.30.31.32 



(7) 



Here ujp is the plasma frequency and has both 

a real and imaginary part related, respectively, to a 
frequency dependent optical effective mass, m*p{u;)/m, 
(where m is the bare electron mass) and scattering rate, 
1/t°p{uj), given by: 



and 



^['m*{uj)/m - f] = -2T,°/{uj) 



(8) 



(9) 



An optical mass renormalization parameter, A°^(w), is 
then defined by m*p{u})/m = 1 + \°p{lo). These quan- 
tities are analogous, but different from, the quasiparti- 
cle mass renormalization and scattering rate which fol- 
low from the quasiparticle self-energy I]''P(a;). Instead, 
Tj°p{uj) is a two-particle quantity^^. It is related to the 
Kubo formula for the current-current correlation func- 
tion which determines cr(w). In general, as stated, the 
optical and quasiparticle mass renormalization and scat- 
tering rate are not the same, although in some limits they 
can be. As an example, the zero frequency mass renor- 
malization = 0) and A^^ are equal. In the end, of 
course, Yi°p{u)) contains exactly the same information as 
does (t(w). In fact, one can solve for 1/t°p and Xop in 
terms of cti and (72 to findi^Si^ 
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47r al(ijj) + cr|(a;) 



and 



1 + A„ 



(72 (w) 



47rw erf (w) -|- ctHuj) 
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(11) 



This does not mean, however, that these quantities have 
no particular value of their own. It is well documented 
that they can speak more directly to certain questions 
than can a itself. For example, 1/t°'p{ijj) provides infor- 
mation on absorption and, in an s-wave superconductor, 
at zero temperature it is zero for photon energies, w, less 
than twice the gap. A coherence peak is seen above this 
energy and at high ld we recover a measure of the nor- 
mal state elastic scattering. Further, it has recently been 
found that in the underdoped cuprates, the real part of 
S°P([j) possesses a "hat like" structure which extends in 
energy over a range of about twice the pseudogap energy. 
This characteristic structure, which is superimposed on a 
large smooth incoherent background is nevertheless un- 
mistakable and provides a direct image of pseuodogap 
formation^ as it enters the in-plane optical data. By 
contrast, the effect of the opening of a pseudogap in the 
real and imaginary part of the conductivity is much more 



subtle and is spread over a much larger energy range. 
Thus, the optical quantities defined in Eq. PH)) and pT|) 
have proved useful and here we will consider how they are 
changed in the case of a distribution of regions of distinct 
gap values. In the top row of Fig. [7] we show results for 
\/T°'P{bj) (solid black curve) in meV as a function of lo 
(also in meV), for the two- and five-gap distributions of 
Figs. |5(b")] and 5(c) The onset of scattering starts at the 
lowest value of twice the gap which is roughly 11 nieV 
in our model. Also shown for comparison (dashed red 
curve) are results for the case of a homogeneous super- 
conductor with the same average gap value. We see that, 
by comparison, the rise in 1/t°p(uj) in the inhomogeneous 
case, is much less steep and proceeds more gradually, be- 
ing spread out over an energy scale which corresponds 
to the spread in gap values in our model distribution. 
The peak in the red dashed curve has its origin in the 
well known coherence peak of the BCS s-wave electronic 
density of states. In the normal state, l/r°P(w) would 
be independent of uj and equal to the input constant im- 
purity scattering rate of 50 meV in this example. This 
is the value to which 1/t°p{uj) (superconducting) tends 
towards for oj greater than a few times the gap. In the 
superconducting state, however, no scattering is possible 
below twice the gap, at which energy it rises sharply and 
overshoots its normal state value, not because the intrin- 
sic scattering potential has increased but rather because 
there are more final states available for scattering. For 
the inhomogeneous case this feature is simply smeared 
out somewhat because of the distribution of different su- 
perconducting regions which are sampled. 

In the lower row of Fig. [7] we show our results for the 
optical mass enhancement parameter [1 -I- \op{uj)\ in the 
two- and five-gap cases. The solid black curve is for the 
inhomogeneous case and the red dashed curve is the ho- 
mogeneous case with an averaged gap value for compar- 
ison. We note the inhomogeneities smear out the peak 
in \op{ijj) around 2A. However, in all other aspects the 
curves are similar. We emphasize two points. First, for u 
large compared with the gap, the mass renormalization 
factor Aop(<j-') tends towards zero, which agrees with its 
normal state value. At w = however, Xop is seen to 
be quite large (~2.5) in both the homogeneous and inho- 
mogeneous cases. In this limit, the optical effective mass 
has a very definite meaning. It is the value of the elec- 
tron mass that one is to give the electrons if one wishes 
to use the classical London formula for the penetration 
depth Ai at T=:0, namely A^^ — Anne^ / m* where n is 
the free electron density, e the electron charge and c the 
velocity of light. In the clean limit, i.e., no elastic scat- 
tering, m* would be the bare electron mass. Including 
impurities in BCS theory reduces the London penetra- 
tion depth and this can be expressed by changing m to 
TO* (a; = 0) in the conventional expression. Our results 
show that this result is only slightly changed in the in- 
homogeneous case, i.e. m*{uj = 0) is only very slightly 
larger. The main effect of inhomogeneities manifests it- 
self at finite frequency in the region of twice the gap. We 
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see that, as compared with the dashed red curve, the sohd 
black curve shows a much broader peak which is a direct 
consequence of spatial gap variations. The peak in X°p{lj) 
around twice the gap is directly tied to the rapid rise in 
1/t°p{uj) at this same frequency, as these are Kramers- 
Kronig (KK) related, i.e., a step in the scattering rate at 
oj = 2A translates into a logarithmic singularity at 2A in 
its KK transform. We note in passing that a recent opti- 
cal study of MoGe (Refi^^) has revealed such a smearing 
of coherence peaks in the thinnest films studied, which 
could be attributed to inhomogeneities, however, other 
evidence suggests otherwise. 
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FIG. 7: (Color online) EMA result for the theoretical self- 
energy of the two grain (left) and five g rain (righ t) composite 
determined by histograms in Figs. 5(b) and 5(c) respectively. 
Top row shows the scattering rate, 1/t°^, and the bottom row 
shows the mass renormalization, 1 -f Aop. 



twice the gap, with the missing optical spectral weight 
transferred to the condensate. For d-wave, the region 
below 2A is also depleted, although some absorption re- 
mains at all frequencies and there is no sharp threshold in 
(Ji{uj) vs to. In this latter circumstance the conductivity 
of the composite system does not exhibit qualitatively 
different behaviour as compared with the single grain 
which would represent a clear signature of the bulk inho- 
mogeneities. In fact, there is no quantitative differences 
between a{Lo) of the composite and the individual grains, 
and the electromagnetic response appears close to that 
for the average gap. By contrast, for an s-wave gap, the 
sharp absorption edge in ai(uj) at twice the gap which is 
characteristic of a uniform BCS superconductor becomes 
less steep reflecting the onset of a distribution of gaps 
rather than of a single one. For a distribution involving 
only two well separated gap values, each is seen distinctly 
in both the real and imaginary parts of the conductivity 
as well as in the real and imaginary parts of the corre- 
sponding optical self energy. This latter property follows 
directly from a knowledge of the conductivity and has 
been found useful in past discussions of, for example, bo- 
son structures in optical properties. For the case of a 
distribution of many gap values, the presence of each in- 
dividual gap is less evident with (Ji{lo) showing rather 
gradual onset over a frequency range representative of 
the variation in gap values. In conclusion, signatures of 
nanoscale spatial inhomogeneities in the far-infrared op- 
tical conductivity are much more readily recognized for 
composites consisting of s-wave gap symmetry compo- 
nents than they are for those with d-wave gap symmetry. 
Indeed, the fact that the d-wave case appears to be un- 
affected by the existence of inhomogeneities in spite of 
their apparent presence in STM, gives support for the 
continued use of single gap models for calculation and 
analysis of optical conductivity in the high-Tc cuprates. 



SUMMARY AND CONCLUSIONS 



Optical absorption is a probe of bulk properties. Using 
effective medium theory, we have calculated the optical 
conductivity of a superconductor consisting of a random 
array of nanoscale regions having different magnitudes 
of superconducting gap. The superconductivity in each 
grain is described within a BCS model, and both s- and d- 
wave symmetry is considered. For simplicity, models con- 
sisting of two and five distinct gap values are examined 
in detail. For an s-wave superconductor, the real part 
of the optical conductivity is zero for frequencies below 
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